Adiabatic and non-adiabatic merging of independent Bose condensates 
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Motivated by a recent experiment [Chikkatur et al. Science, 296, 2193 (2002)] on the merging 
of atomic condensates, we investigate how two independent condensates with random initial phases 
can develop a unique relative phase when we move them together. In the adiabatic limit, the unit- 
ing of independent condensates can be understood from the eigenstate evolution of the governing 
Hamiltonian, which maps degenerate states (corresponding to fragmented condensates) to a single 
state (corresponding to a united condensate) . In the non-adiabatic limit corresponding to the prac- 
tical experimental configurations, we give an explanation on why we can still get a large condensate 
fraction with a unique relative phase. Detailed numerical simulations are then performed for the 
non-adiabatic merging of the condensates, which confirm our explanation and qualitative estima- 
tion. The results may have interesting implications for realizing a continuous atom laser based on 
merging of condensates. 



I. INTRODUCTION 



The prospect of the creation of a continuous, coher- 
ent atomic beam, the atomic analogy of laser beam, 
has been of great interest to many ever since the suc- 
cessful generation of Bose-Einstein condensate (BEC) 
B II HH II II II 0- One of the major challenges to 
such applications lies in the difficulty in continuous con- 
densation of atomic gas due to stringent cooling condi- 
tions. Towards that goal, one needs to spatially separate 
the evaporative cooling from the destructive laser cool- 
ing, which is still pretty challenging even though there've 
been interesting investigations and progresses H, M, US ■ 
Alternatively, one can realize a continuous source of con- 
densate by bringing new condensates into the trap and 
uniting them. Note that independently prepared con- 
densates have a completely random relative phase; in or- 
der to unite them, one needs to have a phase cooling 
mechanism to get rid of this random phase. Interesting 
theoretical investigations have been carried out on the 
possibilities of uniting two existing condensates through 
laser induced phase cooling with condensates confined in 
a high-Q ring cavity pT| . or through phase locking with 
feedbacks from a series of interference measurements |l2| . 

Recently, Chikkatur et al. reported a striking experi- 
ment in which two independently produced condensates 
were merged directly in space by bringing their traps to- 
gether to a complete overlap Jjlj. A condensate fraction 
larger than the initial components was observed from the 
subsequent time-of-flight imaging. This raises the ques- 
tion as how the relative phase between the two compo- 
nent condensates is established during this direct merging 
where there is no additional phase cooling or locking. In 
this paper, we investigate this puzzling phase dynamics 
both in the ideal adiabatic merging limit and in the more 
practical non-adiabatic circumstance. In the ideal adia- 
batic limit, we show that a unique relative phase will be 
established between the two components during direct 
merging, resulting from the eigenstate evolution of the 
governing Hamiltonian (an effective two- mode model). 
The Hamiltonian has degenerate ground states corre- 



sponding to two fragmented condensates when the traps 
are apart. When we overlap the two traps, the degen- 
erate ground states evolve into a single eigenstate corre- 
sponding to a single condensate with a unique relative 
phase between the two initial components. While the 
adiabatic limit requires very slow merging, the real ex- 
periment is actually performed far from that limit |13J . In 
this case, the atoms will slip from the ground state to all 
the nearby eigenstates, and one would not expect to get 
a single condensate. However, through a careful analysis 
of the eigenstate structure of the governing Hamiltonian, 
we argue that one can still get a large condensate frac- 
tion (larger than either of the initial components') with 
a unique relative phase in this non-adiabatic limit. We 
then confirm this evolution picture through detailed nu- 
merical simulations of the merging dynamics. 

We should also point out that the non-adiabatic merg- 
ing considered in this paper, which is based on the pa- 
rameters of the experiment [13 , is not an abrupt connec- 
tion of the two condensates. In the latter case, a unique 
relative phase cannot be established without careful con- 
sideration of dissipation [l^. In the experiment, the time 
scale of the merger( 0.5 s) is long when compared with the 
trap frequency (440Hz) along the merging direction il^. 
So the merging is actually adiabatic with respect to the 
evolution of the condensate wave packets; and dissipation 
due to the quasi-particle excitation in each well should 
not play an important role. However, this same merging 
time scale is very short when compared with the lowest 
excitation energy of the relevant Hamiltonian (the effec- 
tive two-mode model). Many eigenstates of this Hamil- 
tonian will get populated during the merger. That is 
why we call it " non-adiabatic" . Our central observation 
is that due to the particular eigen-energy distribution 
of this relevant two-mode Hamiltonian, even for merg- 
ers much faster than the lowest excitation energy, only 
low-lying states of the Hamiltonian will be populated. 
Furthermore, for all the populated low-lying states, the 
relative phase between the two initial condensate com- 
ponents is well locked. As a result, we will get a large 
merged condensate fraction with a uniquely established 



relative phase. 

In the foUowing, we first present a theoretical model 
for the description of the direct merging process, and 
argue that the internal phase dynamics is mainly cap- 
tured by the simplified two-mode Hamiltonian with time 
dependent parameters, which has been a popular model 
for many theoretical investigations in different scenarios 
[l4, 15, 16, 17, 18, 19]. After that, we go to the main 
topic and investigate the relative phase dynamics. The 
result in the adiabatic limit can be easily understood. In 
the non- adiabatic limit, the arguments need to be based 
on the analysis of the detailed evolution structure of the 
eigenstates of the time-dependent Hamiltonian. Then in 
Sec. Ill, we investigate the relative phase dynamics and 
the evolution of the condensate fraction through detailed 
numerical simulations, with the evolution speed ranging 
from the adiabatic limit to the far non-adiabatic circum- 
stance. We simulate the merging process with different 
configurations and different initial states of the two com- 
ponents. Finally, we discuss the relevance of the results to 
the reported experiment and to the prospect of producing 
continuous atom lasers based on condensate merging. 



II. THEORETICAL MODELING OF 

ADIABATIC AND NON-ADIABATIC MERGING 

OF INDEPENDENT CONDENSATES 

A. An effective two-mode model for the 
description of the relative phase dynamics 

Consider two independently prepared condensates in 
separate traps at low temperature. The two traps are 
then brought closer along a certain direction to allow 
the condensates to merge, as illustrated in Fig. 1. The 
traps can be modeled as a time-dependent double-well 
potential VT(r,t), with the two wells moving towards 
each other. The interactions between the atoms are 
described by the usual (5-pseudopotential, [/(ri,r2) = 
(4:Trhas/m)S{ri,r2), where a^ is the s-wave scattering 
length and m is the mass of the atom. The second quan- 
tization Hamiltonian for such a system has the form 



H 



Ira 



(1) 



where ^I* (r) is the bosonic quantum field operator de- 
scribing the atoms, and the parameter g = ATrhas/m. 

When the traps are not overlapping, the atoms in the 
same condensate are in the same quantum mode. We 
can expand the atomic field operator ^(r) as 'l'(r) = 
ai0i(r) -I- a2(/)2(r),, where 0i(r) (j = 1,2) are the corre- 
sponding mean-fields from the Gross-Pitaevski equation 
for each trap. When we bring the two condensates closer 
by varying the trap potential VT(r,i), we assume that 
the atomic density profile adiabatically follows the move- 
ment of the traps, as is the case in the experiment [13|. 




FIG. 1: Illustration of merging of two independent conden- 
sates. 



With this assumption, at any given time i, we can still 
use a variational mean field 0i(r,t) to describe the cor- 
responding condensate. The condensate wave functions 
now overlap in space, and the atoms in different conden- 
sates can tunnel to each other through the trap barrier. 
As a good approximation, we can expand the atomic field 
#(r,t) at time t as ^(r,t) = ai(/)i(r, t)-|-a202(r,t)- With 
this expansion, the Hamiltonian (1) is simplified to the 
following well-known two-mode model 



fr U{t) f ^ ^ 



J{t) 



{a\a2 + a\ai), (2) 



4 V"^ ' "^ V 2 

where the parameters U (i) and J (t) are given by 



C/(t)-f /(l'/'il' + 



2|0iH02|')rfr, 



J{t) = 



-2J<f>l{^^+V{v,m2dv. 



(3) 



(4) 



In writing Eq. (2), we have assumed a symmetric trap 
with ^ \(t>i\ dr — J 102 1 dr, and we have neglected terms 
that are functions only of the atom number N, as they 
commute with H and can be eliminated by going to a 
rotating frame. We have also neglected the nonlinear 
tunneling terms proportional to {a\a2 + a^ai)'^ (see Ref. 
[l7j for the coefficient before that term) because they are 
either typically smaller than or have the same effect as 
the linear tunneling, so that neglecting them does not 
change the basic physics |23|. Note that given all the 
approximations above, the factor before the last term in 
Eq. (3) has some arbitrariness as this term is on the same 
order of magnitude as some neglected terms in the weak 
overlapping limit. As a convenient choice, we take the 
factor to be 2 so that we have a vanishing U{t) as the 
two traps completely overlap, as one would expect on 
physical grounds. The two- mode Hamiltonian (2), which 
is a very popular model for many theoretical investiga- 
tions in different scenarios [ij, [1^) llfl 113) EM Il3 ' catches 
the main physics of condensates in double-well poten- 
tials resulting from the competition between the tunnel- 
ing and the local nonlinear collision interaction. A more 
detailed and more rigorous derivation of the two-mode 



model from the variational approximation can be found 
in Ref. |l9|. One important limitation for applications 
of the two-mode model is that when the two condensates 
get very close, the coupling between the two condensates 
becomes so strong that some approximations underlying 
the two-mode model are not well justified. In that strong 
overlapping region, Eqs. (2)-(4) serve as a convenient ex- 
trapolation. However, numerical calculation reveals that 
the relative phase and the condensate fraction are es- 
tablished shortly after the commencement of the merger, 
when the two condensates are not yet strongly coupled. 
It is under this context that we think an effective two- 
mode model catches the main physics and is adequate for 
the understanding of the relative phase evolution. 

We will use the time-dependent two-mode model 
(Hamiltonian (2)) to investigate the relative phase dy- 
namics during the merging of independent condensates. 
The temporal behavior of the coefficients U (t) and J (t) 
in principle should be derived from the self-consistent 
variational solutions of (j)i{r,t) and 02(r, t), which is 
pretty involved. Fortunately, in the following we will 
see that the relative phase dynamics is actually insen- 
sitive to the details of U (t) and J (t) . It only depends 
on the rough scale of their rates of variation. Hence, for 
the following discussions in this section, we do not spec- 
ify the explicit forms of U (t) and J{t). Nevertheless, 
we do know the rough trends of the evolutions of U (t) 
and J{t). As the two traps approach each other, J (t) 
will continuously increase from zero to a maximum value 
Jq. Compared with J{t), U (t) typically changes much 
more slowly at the beginning of the merger, and can be 
roughly taken as a constant [/q- As the two traps further 
merge approaching to a complete overlap, one expects 
that U (t) will gradually decrease to zero. These rough 
trends of evolutions of U (t) and J (i) are important for 
the understanding of the relative phase dynamics. 



B. Merging of condensates in the adiabatic limit 

For the convenience of the following discussions, we 
write the two-mode model (2) using effective collective 
spin operators. With the standard Schwinger repre- 
sentation of the spin operators Sx = h{a,\cL2 + aJ^ai), 

Sy — |(a|a2 — 02*^1), and Sz ~ ^{ala2 — ajai), Hamil- 
tonian (2) can be expressed in the following form 



Hit)^U{t)S^^~J{t)Sx. 



(5) 



This Hamiltonian is actually equivalent to that of the 
infinite-range Ising model |2j| , and one has effective anti- 
ferromagnetic interaction as U (t) > 0. 

First, let us assume that the temperature of the atoms 
is low enough, and that the variation speeds of U (t) and 
J (t) are very small. In this ideal adiabatic limit, the 
establishment of a unique relative phase through direct 
merging of condensates can be easily understood from the 
adiabatic theorem, which states that the system in such 



a limit will remain in the ground state of the Hamilto- 
nian H (t). Initially, the two traps are far apart, therefore 
the self interaction among the atoms in each trap domi- 
nates with U (t) -^ Uq and J (t) — > 0. The ground state 
in this case is an eigenstate of S'^ with eigenvalue of 0, 
i.e., with la\a2) = la[ai) — N/2, where N is the total 

atom number. The two condensates have equal number 
of atoms, but there is no phase relation between them, 
which corresponds to fragmented condensates. In the 
final stage of the merger, the coefficients C/ (i) — > and 
J (t) — > Jo J the coupling energy —J (t) Sx thus dominates 
in the Hamiltonian. The ground state in this case is an 
eigenstate of Sx with the largest eigenvalue N/2 (with- 
out loss of generality we assume J (i) > 0). This state 
is actually an eigenstate of the number operator a_^_a-^ of 
the mode a+ — (ai -|- 02) /\/2 with an eigenvalue of N, 
which corresponds to a single condensate. We then have 
a unique zero relative phase between the two initial con- 
densate modes oi and 02, which is due to the fact that 
the coupling energy —J (t) Sx is minimized only when we 
have the same phase between ai and 02. 

The above understanding of the relative phase estab- 
lishment from the adiabatic theorem seems to be simple 
and straightforward, but in a sense it is pretty nontrivial. 
The relative phase really comes from a joint effect of the 
two competitive terms in the Hamiltonian (5). If we have 
only the collision interaction term U (t) Si and start with 
fragmented condensates with (0202) = la\ai\ ~ N/2, 

it's obvious that the condensate fractions will never be 
changed by this collision interaction. On the other hand, 
if we have only the coupling term —J (t) Sx, it will only 
rotate the two condensate modes oi and 02 to two dif- 
ferent modes a'l and a'2, and will never change the con- 
densate fractions either. The condensate fractions can be 
quantitatively described by the eigenvalues of the follow- 
ing single-particle density matrix 19] 



(6) 



The system is in a single-condensate state if p has only 
one macroscopic eigenvalue in the limit N —^ 00; and 
it corresponds to fragmented condensates if p has more 
macroscopic eigenvalues in the large N limit. It is clear 
that the coupling term —J (t) Sx only makes a linear ro- 
tation of the modes ai and a2, and will never change 
the eigenvalues of the single-particle density matrix p. 
Clearly, it is the interplay between these two competing 
interactions in the Hamiltonian (5) that is really respon- 
sible for producing a larger condensate fraction with a 
well-defined relative phase. 

The analysis above gives us a simple understanding 
of how a unique relative phase in principle could be es- 
tablished between two initially independent condensates 
through direct merging in the perfect adiabatic limit. 
However, this analysis hardly explains the merging dy- 




namics in real experiments because the practical config- 
urations are far from this ideal adiabatic limit. Practi- 
cally, first of all, we can never start with an equal number 
state for the two condensates. One can at most control 
the number of atoms in each condensate so that they are 
roughly the same, as each condensate has some number 
fluctuations typically on the order of ^/N . So in reality 
we can not start from the ground state of the Hamiltonian 
(5), but rather a statistical mixture of some low- lying 
states. Secondly, to fulfill the conditions of the adiabatic 
passage, one requires that the variation speeds of U (t) 
and J (i) be significantly smaller than the smallest eigen- 
frequency splitting of the relevant Hamiltonian. The 
Hamiltonian (5) has an energy difference of C/q between 
the ground and the first excited state when U {t) ^ J (t) 
(i.e., when the traps are far apart). This is a tiny energy 
gap, which can be roughly estimated by iJ-c/N, where fic 
denotes the chemical potential of the condensate (~ 2kHz 
|20i) and the atom number is typically around 10^. The 
merging time scale in real experiments (^ 0.5 s 13]) is 
much shorter than the time scale given by I/Uq- This 
shows that we need to go beyond the adiabatic limit to 
understand the merging dynamics in practical configura- 
tions, which is the focus of our next subsection. 



Merging of condensates in the practical 
non-adiabatic circumstance 



If the evolution speed of the Hamiltonian is beyond 
the adiabatic limit, the system prepared in the ground 
state will be excited to some upper eigenstates. The ex- 
citation probability will depend on the ratio of the vari- 
ation rate of the Hamiltonian to the smallest energy dif- 
ference between the two eigenstates J23|- The final state 
in this case is in general a superposition of different low- 
lying eigenstates which have been populated during the 
merger. Therefore, to find out the final state of the sys- 
tem in the non-adiabatic circumstance, it is important to 
look at the evolution of all the eigenstates of the Hamil- 
tonian (5) during the merging process. 

We have calculated the energy spectrum of the Hamil- 
tonian (5) as a function of the ratio of J (i) /U (i) for 
various atom numbers up to a few hundred (see also Ref. 
|25|. where the energy spectrum is drawn for other pur- 
poses). The evolutions of the energy spectrum for dif- 
ferent atom numbers are qualitatively very similar. To 
have a clear display, we only plot in Fig. 2 the whole 
energy spectrum of the Hamiltonian (5) with a small 
atom number N = 20, as the energy levels become too 
dense to be seen clearly in a small figure if N is large. 
Fig. 2 is enough to show some general properties of the 
energy spectrum which arc critical for the understand- 
ing of the merging dynamics outside of the adiabatic 
limit. The spectrum can be divided into three regions, 
corresponding respectively to the Fock, the Josephson, 
and the Rabi region, as discussed in general Josephson 
physics [23. First of all, on the left side of the spec- 



trum (J (i) -^ 0), which corresponds to the Fock region, 
there are a total of iV -I- 1 eigenlevels. While the lowest 
eigenlevel (j = 0) stays non-degenerate, the (2j — l)th 
and (2j)th {j = 1, 2, • • • , N/2) levels are degenerate in 
energy. The two degenerate states can be expressed as 
\N/2+j,N/2-j) and \N/2 - j,N/2 + j) in the Fock 
basis of the modes ai and 02 . The energy difference be- 
tween the ground state and the 2jth eigenlevel is given 
by pU {t) ~ j^Uo- The quadratic scahng of the energy 
difference with j is very important for the following dis- 
cussions. Secondly, on the far right side of the spectrum 
with J ^ NU (not shown in Fig. 2), the system enters 
the Rabi region, where the eigenlevels are almost equally 
spaced. In that region, one can approximately neglect the 
term U{t)J^, and the energy levels are roughly eigen- 
states of the Sx operator. The jth {j = 0, l,2,---iV) 
eigenlevel, expressed as |5 = N/2, Sx — N/2 — j) in the 
Dicke basis of the collective spin operator S , has the 
energy jJ(t) « jJq (compared with the ground state). 
Between these two limiting cases, there is a wide interme- 
diate Josephson region where the two-mode model could 
be understood semi-classically j20j. In the Josephson re- 
gion, with fixed J and [/, the energy levels for the low 
eigenstates are also approximately equally spaced with 
the energy difference between the adjacent levels given 
hy VnUJ. 
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FIG. 2: (Color Online) Energy spectrum of the Hamiltonian 
(5) with the atom number N=20 f25| . 

From the discussions above, we see that the Fock re- 
gion has the smallest energy spacing for the low-lying 
eigenstates with typical experimental parameters. The 
non-adiabatic transition therefore happens dominantly in 
or close to that region. However, due to the quadratic 
scaling of the energy spacing in that region, although 
non-adiabatic transition occurs, the atoms can only be 
excited to some low-lying eigenstates if they start from 
the ground state. The typical merging time in the ex- 
periment T3] is about 0.5 s, significantly longer than the 
time scale of h/ ^c, where fic (^ 2kHz) is the chemical po- 



tential. The parameters NUq and fic are roughly on the 
same order of magnitude for typical atomic gas experi- 
ments. In the Fock region, the energy spacing between 
the 2\/]Vth excited state and the ground state is NUq- So 
if the atoms start from the ground state, the final spread 
in the energy eigenstates due to the non-adiabatic tran- 
sition is at most 2\/N. Of course, in real experiments 
it is not practical to start exactly from the ground state 
\N/2,N/2) of the Hamiltonian (5) when the traps are 
far apart. Let us assume that the initial atom numbers 
Ni and N2 in the two traps are roughly the same with 
Ni ~ N2 — N/2 [231, and that each of them has some 
fluctuations typically on the order of y/N |2J|. Then, 
the initial spread in the energy eigenlevel is about 2\/N. 
Adding together quadratically the contributions from the 
non-adiabatic transfer and from the initial number fluc- 
tuations , we conclude that after the merger, the atoms 
will populate a mixture of the low- lying eigenstates of the 
Hamiltonian (5), with a spread in the eigenstates char- 
acterized by 2y/2N (in terms of the order of magnitude). 
Now we want to show that even if the final state af- 
ter the merger is a mixture of low-lying eigenstates, as 
long as the spread in the eigenstates 2y/2N <C N, we 
still have a large single-condensate fraction. When the 
two traps are merged and we are in the Rabi region, 
we can calculate the condensate fraction in the mode 
Oc — (oi + 02) /V2 for each eigenstate of the Sx oper- 
ator. For the jth eigenstate \S = N/2, S^ = N/2 - j) of 
the Sx operator, its population in the mode Oc is given 
by N — j (This result can be easily seen if we rotate 
Sx to Sz through a linear transformation on the modes 
ai and 02). Now we are in a mixture of the eigenstates 
\S ^ N/2, Sx = N/2 - j) of the Sx operator with j up 
to 2\/2N. For this mixture state, the atom number in 
the particular mode Oc is at least N — 2y/2N. There- 
fore, if the atom number N is huge, we end up with the 
dominant fraction of atoms in the mode Oc even if the 
merging is far outside the perfect adiabatic limit. The 
basic reason for this result lies in the fast growth of the 
excitation energy with the levels. The system is kept in 
some low-lying states of the Hamiltonian (5) although 
the evolution is non-adiabatic. When most of the atoms 
are in the condensate mode Oc, a unique relative phase 
of zero is correspondingly established between the two 
initial modes oi and 02. 



we can almost get a single condensate even if the merg- 
ing is far outside the adiabatic limit {t,n ^ h/Uo). In 
this section, we would like to quantitatively test this re- 
sult through more detailed numerical simulations of the 
merging dynamics. 

To quantitatively describe the merging dynamics, we 
look at the evolution of the largest condensate frac- 
tion derived from the single-particle density matrix (6). 
The largest eigenvalue of the single-particle density ma- 
trix p can always be expressed as an expectation value 
(^aj.ac), where Oc is a rotated mode corresponding to 
the largest condensate fraction and is defined as Oc — 
cos9ai + sin0e*'^a2. The rotation angle 6 reflects the 
contributions to the new condensate mode Oc from the 
two initial modes ai and 02, and ip specifies their relative 
phase in this new mode. With this characterization of 
the relative phase, here we ignore explicit discussion on 
the inherent quantum uncertainty of the phase operator 
caused by the finite total atom number. This quantum 
uncertainty is on the order of 1/\N , where the atom 
number N ^ 10^ for typical experiments, so it repre- 
sents a small effect. We refer the readers to Ref. |2g 
for detailed discussions on that issue. Numerically, we 
start with certain initial states of the modes oi and 02, 
as will be detailed in the following subsections, solve the 
evolution of these states under the Hamiltonian (5) , and 
calculate the single-particle density matrix for each in- 
stantaneous state to find out the evolution of the atom 
number {a\ac) in the largest condensate fraction mode 
Qc and the corresponding parameters 6 and Lp. The atom 
number (ajoc) is given by the largest eigenvalue of p, 
while 9 and (p are specified by the corresponding eigen- 
vector (cos 0, sin 06*'''). 

To quantify the state evolution, we need to specify the 
parameters U (t) and J (t) in the Hamiltonian (5) as a 
function of time. To be consistent with the experimental 
configurations in Ref. (13), we assume that the two initial 
condensates are confined in cigar-shaped traps which are 
moving towards each other along the radial direction with 
uniform speed. As we have mentioned before, the details 
of the functions U (t) and J (t) are not important, but 
rather their time scales. So, for simplicity, we assume 
that the wave functions 0i(r,t) {i — 1,2) of the modes 
Qi , which adiabatically follow the movement of the cigar- 
shaped traps, also have cigar-shaped profiles specified by 
the following Gaussian function 



III. NUMERICAL SIMULATION OF THE 

MERGING DYNAMICS OF TWO 

CONDENSATES 

A. Numerical simulation methods 

In the previous sections, we have shown that through 
a direct merger it is possible to establish a unique rel- 
ative phase between two initially independent conden- 
sates, and argued that this relative phase is robust and 
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where Up and Cz are the widths of the wave functions 
along the radial and the axial directions, respectively. 
The centers of the wavepackets are given by a;io(i) = 
xo (1 - t/t„i) /2 and X2o{t) = -Xq (1 - t/tm) /2, defining 
a constantly diminishing separation from xq to through- 



out the merging time tm- The variation of the parame- 
ters U (i) and J {t) are then calculated with Eqs. (3)(4) 
from these postulated wave functions. The typical evo- 
lution of U (t) and J (t) are shown in Fig. 5(a-b) with 
az = lOiTp, following the experimental conditions [l3J. 
Although in real experiments the total atom number is 
normally on the order of '^ 10^, in numerical simulations, 
due to the limited computation speed, we can only deal 
with moderate atom numbers with A'' ~ 10^. However, 
we know that NUq and Jq are typically on the same order 
of magnitude with Jq somewhat smaller than NUq. In 
order to be consistent with the practical configurations, 
we should re-scale the parameters U (t) and J (t) with a 
global constant so that NUq ^ Jq. In our simulation, 
without loss of generality, we assume NUq = AJq. When 
U (t) and J (i) are specified, we can numerically calculate 
the evolution of the atom number {a].a^ in the largest 
condensate fraction and the relative phase ip with the 
method described above. 



B. Merging of two condensates in Fock states 

We first investigate the merging of two condensates 
in Fock states, with the atom numbers A^i = 51 and 
N2 = 49, respectively. The numerical results are shown 
in Figs. (3a)- (3h), where the largest condensate fraction 
77 (defined as 77 = (ajoc) /N with N = Ni + N2) and 
the parameters {O^Lp) are plotted against the dimension- 
less time UqI. Note that for the initial Fock state, the 
relative phase (/? is not well defined (random). In the 
numerical program, one basically randomly picks up a 
particular initial phase. In Fig. 3, the initial value of Lp 
is set to zero simply due to the convention of the pro- 
gram (it sets phase (,c = for the complex number Oe*"^). 
The merging dynamics and the final relative phase do 
not depend on this particular choice. We test other as- 
signments of the initial phase, and basically there are 
no changes in the dynamics (except for different initial 
jumps of the phase which are meaningless). We have 
different merging time scales for figures (a)-(h), ranging 
from the adiabatic limit C/oim = 40 to the complete non- 
adiabatic limit f/o^m = 0.04, where tm denotes the time 
of the merger. With Uotm = 40 :^ 1, the picture of the 
adiabatic mapping is basically valid, and as expected, 
finally we have more than 95% of atoms evolving into 
the largest condensate fraction (see Fig. 3g). For Figs. 
(3e) and (3c), the time scales change from t/oim = 4 to 
Uoitm = 0.4, which is certainly outside the adiabatic limit 
as f/o^m £ 1- However, since NUotm = lOOUotm ^ 1, as 
we analyzed in Sec. IIC, only some low-lying eigenlevels 
will be populated during the evolution, and we should 
have a good condensate fraction. This is confirmed by 
the figures (3e) and (3c). Compared with Fig. (3g), the 
final condensate fraction is reduced a little bit, but not 
much, in particular for Fig. (3e) with NUotm 3> 1. Fi- 
nally, in Fig. (3a) the variation is so fast {NUotm — 4) 
that many of the eigenlevels of the Hamiltonian (5) will 



be populated and we expect to have a poor final conden- 
sate fraction. This is confirmed by the numerical simula- 
tion, where we can hardly see any considerable increase 
of the condensate fraction. 

Figures. (3b) (3d) (3f) and (3h) show the evolution of 
the corresponding condensate mode which has the largest 
fraction of the atoms. For Figs. (3d) (3f) and (3h), we 
expect that the final state should be close to the low- 
est eigenstate of the operator ~Sx- The corresponding 
condensate mode should be Qc — (ai -I- 02) /\/2, with 
9 — > 7r/4 and 1^9 —> 0. Indeed, we can see from these three 
figures that the parameters 9 and ip approach their sta- 
tionary values with very small oscillations as soon as J (t) 
becomes comparable to U{t). When J (t) ~ U{t), the 
system is well inside the Josephson region. This means 
that a unique relative phase has been established between 
the two initial modes ai and 02 when their coupling is still 
weak. In such a region, the two-mode model is still well 
justified, which supports the application of the two-mode 
model in understanding the establishment of the relative 
phase. The signature from Fig. (3b) is less clear because 
the corresponding increase in the condensate fraction still 
remains to be small. 



C. Merging of t^vo condensates in Fock and 
coherent states 



In this subsection, we investigate the merging of two 
condensates in a Fock state with TVi = 50 and in a co- 
herent state |a) with |a| = 64, respectively. This is 
motivated, on the one hand, by the consideration that 
the source for an operating atom laser resembles a co- 
herent state; and on the other hand, by the curiosity to 
find out the influence of the initial number fluctuations 
on the merging dynamics. A coherent state is certainly 
not an eigenstate of the Hamiltonian (5) at i = with 
J (t) — > 0, so we start from a superposition of a series of 
eigenlevels instead of a particular one. 

The simulation results are shown in Figs. (4a)- (4e) for 
the evolutions of the largest condensate fraction rj and 
the corresponding parameters {9, p), with the merging 
time varying from Uotm = 0.04 to Uotm = 4. The results 
are qualitatively similar to the ones displayed in Fig. 3. 
The main difference is that the final condensate fraction 
corresponding to the same time scale is notably worse 
now. This is understandable, as the initial state we start 
with is not an eigenstate of the governing Hamiltonian 
and it has significant initial number fluctuation. The co- 
herent state has number fluctuation on the order of vN, 
so from the analysis in Sec. IIC, we expect that the final 
condensate fraction may reduce by an amount of order 
y/N/N, which is pretty close to the results shown in Figs. 
(4a)-(4e). In our simulation, the coherent state is not 
large due to the limitation of the computation efficiency; 
1/\/N is therefore not a small factor and we see con- 
siderable decrease in the condensate fraction. In more 
practical cases with N ~ 10®, the initial number flue- 
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FIG. 3: (Color Online) Evolution of the largest condensate fraction and the relative phase parameters for various merging 
time scales represented by the dimensionless parameter Uotm- Figs. (a),(c),(e),(g) are for evolution of the largest condensate 
fractions, with the merging time scales Uotm = 0.04, Uotm ~ 0.4, Uotm ~ 4, Uotm = 40, respectively. Figures (b),(d),(f),(h) 
are for evolution of the rotation angle 6 (solid line) and the relative phase ip (dashed line), with the merging time scales 
corresponding to Figs. (a),(c),(e),(g), respectively. 



tuation of each condensate should have a much smaller 
influence on the final condensate fraction. 



D. Comparison of different merging methods 

We have mentioned before that the merging dynamics 
is only determined by the rough time scale in the vari- 
ations of U (t) and J (t) , rather than the explicit func- 



0.58 



0.56 



O 


1 










^ 


OQ 


^ 




Hh 




D 


OX 






K! 




^ 




C 


07 


o 




T3 




o 


0.6 


U 






0.5 




1 




0.9 




0.8 




0.7 




0.6 




0.5 




> 



Pi 



0.2 0.3 0.4 




2 
1.5 

1 
0.5 











■\ ■ 


1 \ 


^^^-^^^ 




1 \ 




~ ^-^^^^ 


N, 




' 


1 


~ ^ — _ _ _ 


^ — — — " 


1 

1 




(b) 







0.01 0.02 0.03 0.04 




m 



m 



FIG. 4: (Color Online) Similar to Fig. 3, but with different states for the initial component condensates. One of them is in a 
Fock state with N=50 and the other is in a coherent state \a) with a = 8. The merging time scales are given respectively by 
UoU, = 0.04 (Figs. (a),(b)), (7ot™. = 0.4 (Figs. (c),(d)), C/oU = 4 (Figs. (e),(f)) 



tional forms of these two parameters. We now verify this 
result by considering different merging methods. Let us 
consider two cigar-shaped condensates (with a^ — lOcTp 
in Eq. (7)) merged along the radial and the axial direc- 
tions, respectively. The variations of U (t) and J (t) as 
functions of time are shown in Fig. (5a) and (5b) for 
these two cases. For merging along the axial direction, 
the evolutions of the largest condensate fraction -q and the 
corresponding mode parameters {9,(p) are shown in Fig. 
(5c) and (5d) with Uatm = 4. Compared with the corre- 
sponding results for merging along the radial direction, 
we see there is little difference in the final condensate 
fraction although the variations of U (t) and J (t) dif- 
fer considerably in the two cases. This shows that what 



matters most for the merging dynamics is the rough time 
scale of U (t) and J {t) as we have mentioned before. It 
also justifies our use of the simple Gaussian functions in 
Eq. (7) to model the individual condensate wave func- 
tions, as the functional details of U (t) and J {t) are not 
so important. 

This result should not be misunderstood as that it is 
equally good in real experiments to merge the conden- 
sates along the radial or the axial direction. Experimen- 
tally, for the cigar-shaped condensates, it is much better 
to merge them along the radial direction |l3|. The rea- 
son is that we have assumed that the density profile of 
the atoms can adiabatically follow the movements of the 
traps to validate the two-mode model. Although that is 




J(t) 0.05 




Tl 



0.8 



0.7 



0.6 



0.5 





p/r^""^ 


' 


^ _ ^ ^ ^ 






'1 


1 ^ 






j 






(c) 



1 2 

m 






Pi 




FIG. 5; (Color Online) Comparison of the dynamics for merging along different directions. (a)(b) U and J as functions of the 
dimensionless time while merging along the radial direction (solid line) or the axial direction (dashed line).(c) Comparison of 
the evolution of the largest condensate fraction r). The solid line represents merging along the radial direction and the dashed 
line represents merging along the axial direction. Despite different variations of U and J shown in (a)(b), there seems to be 
little difference in the evolution of the condensate fraction, (d) Evolution of the rotation angle ^(solid line) and the relative 
phase (/^(dashed line) when merging along the the axial direction. Qualitatively, it is similar to the evolution shown in Fig. 3(f) 
for merging along the radial direction. 



the case for merging along the radial direction, it will be 
much harder to meet this condition if the condensates are 
merged along the axial direction. In the latter case, one 
needs to have the condensates further away from each 
other initially to have negligible J {t) at the beginning; 
one also needs to move the atoms significantly faster to 
have the same time scale for merging. However, the trap 
along the axial direction is much weaker, and a large frac- 
tion of the atoms could be excited during the merger. As 
a result, the atomic density profile would be left behind 
the movements of the traps unless one reduces the merg- 
ing speed to an undesirable value (for instance, with a 
time longer than the condensate life time). Therefore, 
merging along the weaker trapping direction makes the 
individual condensate wave functions harder to follow 
adiabatically the movements of the traps, which would 
invalidate the two-mode approximation. Nevertheless, if 
the two-mode approximation could be justified, the rel- 
ative phase dynamics would then become insensitive to 
the detailed merging methods, as shown by this numeri- 
cal simulation. 



IV. SUMMARY 



Using a two-mode model, we have studied the dynam- 
ics of relative phase and condensate fraction during the 
direct merging of two independently prepared conden- 
sates with random relative phase. In accordance with 
a recent experiment [13|, we found that it is possible to 
create a single condensate with larger condensate frac- 
tion and a unique zero relative phase between the initial 
condensate modes. By examining the energy spectrum of 
the Hamiltonian under the two-mode approximation, the 
process can be understood both within and without the 
adiabatic limit. In the ideal adiabatic limit, the result 
can be explained using the adiabatic theorem, and the 
system will evolve from a fragmented condensate to a sin- 
gle condensate following the evolution of the eigenstate 
of the Hamiltonian during the merger. Beyond the adia- 
batic limit, careful analysis of the evolution of the eigen- 
spectrum is needed. Qualitatively, due to the quadratic 
scaling of the excitation energy with the energy levels, 
only the low-lying eigenstates of the Hamiltonian will be 
populated even if the process is far from the adiabatic 
limit. At the end of the merger, the mixture of those 
low-lying states will give rise to the final state with dom- 
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inant number of atoms in the desired single condensate 
mode. Numerical simulations are then performed, and 
the results are in good agreement with our analysis. The 
results may have interesting implications for realization 
of a continuous atom laser based on direct merging of in- 
dependently prepared condensates. Because of the simi- 
larity of our model to the case of condensates in double 



wells or optical lattice, the adiabatic and non-adiabatic 
evolution picture described here may also find applica- 
tions in controlling the dynamics of condensates in those 
optical potentials. 
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